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ABSTRACT 

During  the  Higher  Harmonic  Control  Aeroacoustic  Ro¬ 
tor  Test,  extensive  measurements  of  the  rotor  aerody¬ 
namics,  the  far-field  acoustics,  the  wake  geometry  and 
the  blade  motion  for  powered,  descent,  flight  condi¬ 
tions  were  made.  These  measurements  have  been  used 
to  validate  and  improve  the  prediction  of  blade- vortex 
interaction  (BVI)  noise.  The  improvements  made  to 
the  BVI  modeling  after  the  evaluation  of  the  test  data 
are  discussed.  The  effects  of  these  improvements  on 
the  acoustic-pressure  predictions  are  shown.  These 
improvements  include  re-structuring  the  wake,  mod¬ 
ifying  the  core  size,  incorporating  the  measured  blade 
motion  into  the  calculations  and  attempting  to  im¬ 
prove  the  d3nnamic  blade  response.  A  comparison  of 
four  different  implementations  of  the  Ffowcs  Williams 
and  Hawkings  equation  is  presented.  A  common  set  of 
aerodynamic  input  has  been  used  for  this  comparison. 

1  INTRODUCTION 

In  order  to  correctly  predict  the  aerodynamic  and 
acoustic  pressure  for  a  helicopter  operating  in  pow¬ 
ered,  descent  flight,  it  is  of  utmost  importance  to  cor¬ 
rectly  model  the  wake  geometry  of  the  vortices  shed 
from  the  rotor  tips.  The  wake-geometry  model  should 
contain  information  on  the  vortex  location,  the  vor¬ 
tex  strength  and  the  vortex  core  size.  In  a  recent, 
international  cooperative  test.  Higher  Harmonic  Con¬ 
trol  (HHC)  Aeroacoustic  Rotor  Test  (HART),  exten- 
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sive  measurements  of  the  rotor  aerodynamics,  the  far- 
field  acoustics,  the  wake  geometry,  and  the  blade  mo¬ 
tions  were  made.^  In  preparation  for  this  test,  an 
international  team  of  researchers  from  DLR,  Germany, 
ONERA,  France,  the  US  Army  Aeroflightdynamics  Di¬ 
rectorate  (AFDD),  and  NASA  Langley  have  been  im¬ 
proving  the  prediction  capability  of  blade  vortex  in¬ 
teraction  (BVI)  phenomena.  This  prediction  capabil¬ 
ity  includes  models  of  the  vortex  geometry,  aeroelastic 
blade  motion,  aerodynamic  surface  pressure,  and  far- 
field  acoustic  predictions. 

The  first  validation  effort  of  the  prediction  team  was 
documented  in  Ref.  2.  This  validation  effort  included 
wake  geometry  predictions,  aerodynamic  surface  pres¬ 
sure  predictions  and  acoustic  far-field  predictions  of 
the  2-bladed,  1/7  scale  model  OLS  rotor  tested  in 
the  Duits-Nederlandse  Windtunnel  (DNW)  in  1982.^*^ 
Leading  edge  pressure  transducers  at  x/c  =  0.03  and 
far-field  acoustic  data  from  this  test  were  used  for  the 
validation  efforts.  The  prediction  codes  compared  rea¬ 
sonably  well  with  the  test  data  and  captured  the  qual¬ 
itative  characteristics  of  the  acoustic  test  data.  How¬ 
ever,  further  improvements  to  the  wake  geometry  mod¬ 
els  are  necessary  to  adequately  predict  the  quantita¬ 
tive  characteristics  of  the  acoustic  pressure  needed  to 
evaluate  noise  reduction  techniques. 

The  second  validation  effort  of  the  HART  prediction 
team  focused  on  the  prediction  of  the  wake  geom¬ 
etry,  aeroelastic  blade  motion,  aerodynamic  surface 
pressure,  and  far-field  acoustics  for  the  dynamically 
scaled,  four-bladed,  40%  model  of  the  bearingless  BO- 
105  main  rotor.^  The  predictions  are  documented  in 
Ref.  6.  Unfortimately,  only  the  far-field  acoustic  pres¬ 
sures  were  available  for  this  validation  effort.  There¬ 
fore,  the  wake  geometry,  the  aeroelastic  and  the  aero- 
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dynamic  predictions  of  each  prediction  team  member 
were  compared  to  one  another  to  provide  a  better  un¬ 
derstanding  of  how  each  method  predicted  the  BVI 
phenomena.  This  exercise  proved  the  importance  of 
including  the  aeroelastic  effects  in  the  prediction  of 
BVI  noise.  However,  the  most  important  conclusion 
of  this  exercise  was  that  the  prediction  team  needed 
to  improve  the  modeling  of  the  effect  of  4/rev  HHC  on 
BVI  noise  for  the  BO105  model  rotor. 

The  acoustic  predictions  of  the  third  validation  effort 
are  documented  in  the  current  paper.  The  third  val¬ 
idation  effort  is  focused  on  predicting  test  points  run 
during  the  HART  program.  The  blades  used  for  this 
test  are  the  same  BO105  model  scale  blades  used  for 
the  European  HELINOISE  aeroacoustic  rotor  test.*^ 
Both  the  pre-test  and  the  post-test  acoustic  predic¬ 
tions  are  presented.  The  pre-HART  acoustic  predic¬ 
tions  are  shown  to  indicate  the  status  of  the  prediction 
ability  before  the  HART  data  was  available  to  help  the 
prediction  team  improve  their  BVI  phenomena  model¬ 
ing  capability.  The  pre-HART  predictions  were  done 
to  help  define  the  test  matrix,  improve  the  measure¬ 
ment  techniques  and  to  provide  background  for  possi¬ 
ble  unexpected  experimental  results.  A  brief  summary 
of  the  improvements  made  to  the  BVI  modeling  will 
also  be  included.  For  more  information  on  the  details 
of  the  improvements  to  the  BVI  modeling,  see  Ref.  8. 

All  the  prediction  team  members  use  a  numerical  sim¬ 
ulation  of  the  Ffowcs  Williams  and  Hawkings  (FW-H) 
equation^  to  predict  the  acoustics.  Therefore,  all  pre¬ 
dictions  require  detailed  information  about  the  aero¬ 
dynamic  surface  pressure  on  the  rotor  blade.  The  aero¬ 
dynamic  surface  pressure  calculations  for  the  acoustic 
results  shown  in  this  paper  are  documented  in  Ref.  8. 
In  order  to  directly  compare  the  different  implemen¬ 
tations  of  the  FW-H  equation  of  each  team  member, 
a  common  set  of  blade  surface  pressures  was  used  by 
all  team  members.  These  comparisons  document  any 
differences  in  acoustic  predictions  due  to  the  different 
implementations  of  the  FW-H  equation. 

2  THE  ACOUSTIC  PREDICTION  CODES 

For  the  acoustic  calculations,  the  DLR  has  developed 
the  acoustic  code  AKUROT,^®  ONERA  computes  the 
acoustics  using  the  PARIS  code^^,  AFDD  utilizes  the 
RAPP  codei2  p^aSA  Langley  runs  WOPWOP^^ 
to  predict  the  acoustics.  All  the  acoustic  prediction 
codes  use  only  the  linear  thickness  and  loading  terms  of 
the  Ffowcs  Williams  and  Hawkings  equation,  assuming 
that  the  noise  from  the  nonlinear  quadrupole  term  is 
negligible  for  the  blade- vortex  interactions  studied  in 
this  program. 

2.1  DLR  Acoustic  Prediction  Methodology 

To  predict  the  acoustic  pressure,  the  DLR  uses  the 
acoustic  code  AKUROT.  This  code  was  developed  in 
1987^  and  first  validated  with  experimental  data  from 


the  1982  AHl-OLS  model  rotor  test  in  the  DNW.  The 
code  is  based  on  Farassats  formulation  1  derived  from 
the  FW-H  equation  and  includes  a  thickness  noise  term 
and  a  loading  noise  term  for  the  HART  calculations. 

The  input  data  for  the  loading  term  are  provided  by 
the  S4-rotor  simulation  code.®  Since  the  aerodynamic 
model  in  the  S4  code  is  based  on  lifting  line  theory, 
only  sectional  loading  data  is  available  at  20  radial 
stations  and  180  azimuthal  stations.  For  \ise  in  the 
acoustic  prediction  code,  the  loading  data  were  inter¬ 
polated  to  get  inputs  for  1024  azimuthal  and  160  radial 
stations.  The  acoustic  code  uses  a  chordwise  compact 
model  of  the  loading  noise. 

2.2  NASA  LaRC  Acoustic  Prediction 
Methodology 

NASA  Langley  uses  the  rotor  acoustic  code 
WOPWOP^^  which  implements  the  acoustic  for¬ 
mulation  lA  of  Farassat  to  predict  the  acoustic 
pressure.  The  input  to  WOPWOP  includes  the 
physical  characteristic  of  the  rotor  blade  and  the 
aerodynamic  blade  loading  as  a  function  of  az¬ 
imuthal  position.  For  the  predictions  presented, 
the  loads  were  input  to  WOPWOP  every  1  de¬ 
gree  of  azimuth.  When  the  aerodynamic  input 
used  in  WOPWOP  is  calculated  by  HIRES^^,  the 
loading  noise  is  modeled  by  a  chordwise  compact 
source.  When  a  chordwise  distribution  of  the  blade 
loading  is  available  from  the  fiill  potential  code 
FPRBVI,  the  loading  noise  is  modeled  by  noncompact 
sources  distributed  over  the  chord. 

2.3  ONERA  Acoustic  Prediction  Methodology 

At  ONERA,  the  noise  radiation  is  computed  by  the 
PARIS  code^^  using  a  pressure  distribution  provided 
by  ARHIS.^^  The  azimuthal  step  size  of  the  input  data 
depends  on  the  impulsivity  of  the  interactions.  An  ef¬ 
ficient  spanwise  interpolation  model  has  been  imple¬ 
mented  in  PARIS  in  order  to  minimize  the  amount 
of  airload  data  required  for  BVI  noise  predictions. 
This  interpolation  method  identifies  the  BVI  impul¬ 
sive  events  in  the  acoustic  signatures  generated  at  each 
individual  blade  section  and  accounts  for  their  phase 
and  amplitude.  This  method,  used  with  aerod3mamic 
data  from  10  blade  sections,  provides  the  same  result 
as  a  prediciton  without  interpolation  that  uses  airload 
data  from  30  blade  sections.  PARIS  provides  the  cor¬ 
relation  between  the  peaks  of  a  BVI  noise  signature 
and  the  corresponding  source  locations  on  the  rotor 
disk. 

2.4  AFDD  Acoustic  Prediction  Methodology 

At  AFDD,  RAPP^^,  a  rotor  acoustic  prediction  pro¬ 
gram  based  on  the  FW-H  equation,  calculates  the 
acoustic  pressure  using  the  FPR^®”^®  computed  sur¬ 
face  pressures.  RAPP  can  use  either  the  chordwise 
compact  or  noncompact  formulation.  When  the  force 
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terms  in  the  FW-H  equation  are  modeled  as  chordwise 
compact  sources,  the  sources  are  distributed  spanwise 
along  the  quarter  chord  of  the  acoustic  planform.  The 
acoustic  planform  consists  of  the  locations  of  the  con¬ 
tributing  sources.  This  method  is  referred  to  as  the 
acoustic  lifting  line  method. 

When  the  noncompact  formulation  is  used  to  model 
the  force  terms,  20  sources  are  distributed  over  the 
airfoil  surface  at  each  spanwise  location.  For  both  for¬ 
mulations,  11  spanwise  stations,  distributed  over  the 
outer  60%  of  the  span,  are  used  in  the  computation. 
The  time  step  for  the  acoustic  calculation  is  0.7  de¬ 
grees.  The  time  resolution  of  the  input  aerodynamic 
loading  is  0.125  degrees.  The  advantage  of  the  com¬ 
pact  formulation  as  compared  to  the  noncompact  for¬ 
mulation  is  the  reduced  time  and  memory  required. 
For  the  compact  formulation,  only  the  section  loading 
is  needed  at  each  spanwise  location  for  each  time  step 
(0.125  degrees).  Whereas,  for  the  noncompact  formu¬ 
lation,  the  pressure  at  the  20  surface  points  is  needed 
at  each  spanwise  station  for  each  time  step. 

3  TEST  AND  MODEL  DESCRIPTION 

The  HHC  aeroacoustic  rotor  test  configuration  con¬ 
sisted  of  a  Mach  scaled,  40%  model  of  the  hingeless 
BO  105  main  rotor  mounted  on  the  DLR  test  rig.  The 
rectangular  rotor  blades  have  -8  degrees  of  linear  twist, 
a  radius  of  2  meters,  a  constant  chord  of  0.121  me¬ 
ters  and  a  modified  NACA  23012  airfoil.  The  model  is 
highly  instrumented  with  124  pressure  transducers  and 
53  strain  gages.  The  test  was  run  in  the  6m  x  8m  open 
jet  configuration  of  the  DNW.  The  acoustic  measure¬ 
ments  were  made  with  an  array  of  eleven  microphones 
mounted  on  a  ground  based  traverse  mechanism.  The 
microphones  are  arranged  symmetrically  with  respect 
to  the  rotor  hub  and  are  spaced  0.54  meters  apart.  The 
traverse,  located  2.3  meters  below  the  hub  of  the  rotor, 
travels  from  4  meters  upstream  of  the  rotor  hub  to  4 
meters  downstream  of  the  rotor  hub.  Data  is  acquired 
while  the  microphone  array  is  continuously  traversing 
downstream,  but  the  actual  measurement  locations  are 
at  0.5  meter  intervals  starting  at  4.0  meters  upstream 
of  the  hub  resulting  in  187  measurement  locations  be¬ 
low  the  rotor  hub.  Figure  1  shows  the  resulting  grid  of 
measurement  locations  with  respect  to  the  rotor  plane. 
One  microphone  was  located  on  the  body  of  the  test  rig 
and  two  microphones  were  mounted  in  the  plane  of  the 
rotor  on  the  8m  x  6m  nozzle.  For  each  test  condition, 
the  acoustic  pressure  was  recorded  for  30  rotor  revolu¬ 
tions  at  a  rate  of  2048  data  points  per  revolution.  The 
pressure  time  histories  and  sound  pressure  level  con¬ 
tours  presented  in  this  paper  for  code  validation  are 
based  on  an  ensemble  average  of  30  revolutions. 

The  HART  test  matrix  consisted  of  descent  flight  con¬ 
ditions  with  a  variation  of  the  glide  path  angle  from  -3 
degrees  to  -1-12  degrees  and  a  variation  of  the  advance 
ratio  from  0.144  to  0.275.  The  nominal  test  condition, 
or  the  baseline  case,  simulated  a  landing  approach  with 


2.7  4 


Fig.  1  Grid  of  187  microphone  locations  below  the 
rotor  plane. 

a  glide  path  angle  of  6  degrees  and  an  advance  ratio, 
/i,  of  0.15.  This  case  corresponds  to  a  shaft  tilt  angle, 
a„  of  5.3  degrees  and  a  wind  tunnel  velocity  of  33  me¬ 
ters/second.  For  this  baseline  condition,  a  systematic 
variation  in  phase  shift  of  the  HHC  was  performed  for 
3- ,4-  and  5/rev  input  of  0.8  degree  amplitude.  The 
most  effective  HHC  input  for  reducing  noise  and  vi¬ 
brations  was  the  3/rev.  Therefore,  the  measured  data 
and  the  analytic  solutions  presented  in  this  paper  for 
code  validation  will  be  the  baseline  case  and  2  cases 
with  3/rev  HHC  input  of  0.8  degrees.  One  case  corre¬ 
sponds  to  low  noise  and  one  case  corresponds  to  low 
vibrations.  The  test  conditions  presented  in  this  paper 
are  listed  in  Table  1. 

4  RESULTS 

For  the  test  cases  listed  in  Table  1,  the  measured 
and  predicted  acoustic  pressure  time  histories  are  pre¬ 
sented  at  the  location  of  maximum,  mid-frequency, 
sound  pressure  level  (SPL)  on  the  retreating  side  and 
on  the  advancing  side  of  the  rotor.  The  locations  of 
maximum,  mid-frequency,  SPL  are  different  for  each 
test  condition  and  are  listed  in  Table  2.  The  average 
time  history  at  each  measurement  point  below  the  ro¬ 
tor  (see  Fig.  1)  is  used  to  determine  the  mid-frequency 
SPL.  The  band  limit  is  the  6th  to  the  40th  blade  pas¬ 
sage  frequency.  This  mid-frequency  band  is  believed 
to  contain  the  frequencies  present  in  BVI  noise.  The 
measured,  mid-frequency,  SPL  contour  plots  for  the 
three  test  cases  listed  in  Table  1  are  shown  in  Fig. 
2.  The  measured,  acoiistic-pressure  time  histories  are 
shown  in  Fig.  3  at  the  maximum  SPL  locations  listed 
in  Table  2.  Please  note  that  in  normalized  time,  one 
rotor  revolution  is  equal  to  one  unit  of  time. 

The  pre-test  and  post^test  acoustic  predictions 
(acoustic-pressure  time  histories  and  mid-frequency, 
SPL  contours)  from  each  research  center  are  presented 
in  the  following  sections  to  indicate  which  improve¬ 
ments  to  the  BVI  prediction  methodology  have  the 
greatest  potential.  A  brief  description  of  the  changes 
made  to  the  BVI  prediction  methodology  by  each  re- 
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x(m) 


Cj’ 

a,  (deg) 

RPM 

HHC  amp  (deg) 

HHC  phase  (deg) 

Baseline  (Rim  140) 

0.00432 

0.152 

5.3 

1041 

0.0 

0.0 

Low  Noise  (Run  138) 

0.00413 

0.152 

5.3 

1041 

296 

Low  Vibration  (Run  133) 

0.0044 

0.151 

5.3 

1041 

177 

Table  1  Test  points  used  for  code  validation 


advancing  side 

retreating  side 

X  (m) 

y  (m) 

z  (m) 

X  (m) 

y  (m) 

z  (m) 

Baseline  (Run  140) 

0.0 

2.18 

-2.29 

2.0 

-1.07 

-2.28 

Low  Noise  (Run  138) 

-0.5 

1.1 

-2.29 

2.0 

-1.07 

-2.28 

Low  Vibration  (Rim  133) 

0.05 

2.18 

-2.29 

2.5 

-1.07 

-2.28 

Table  2  Microphone  positions  for  maximum  SPL 


Fig.  2  Measured,  mid- frequency  (6th-40th  harmonic),  SPL  contour  plots  for  the  test  cases  listed  in  Table  1 


Normalized  time 


Low  Noise  (Run  138)  Low  Vibration  (Run  133) 


Normalized  time  Normalized  time 


Fig.  3  Measured,  acoustic-pressure  time  histories  for  the  test  cases  listed  in  Table  1.  Row  1)  advancing  side 
microphone  locations  listed  in  Table  2.  Row  2)  retreating  side  microphone  locations  listed  in  Table  2. 


search  center  is  included. 

4.1  DLR  Pre-HART  and  Post-HART  Acoustic 
Predictions 

The  pre-HART  acoustic  predictions  made  by  the  DLR 
are  based  on  spanwise  blade  loadings  provided  by  the 
S4  rotor  simulation  code  described  in  Ref.  6.  For  the 
post^HART  acoustic  predictions,  the  wake  model  in 
S4  has  been  modified.  The  modification  to  the  wake 
model  is  mainly  an  application  of  a  wake  contraction 
model.  A  15%  contraction  has  been  applied  to  the  to 
the  baseline  case  and  low  noise  case  and  a  10%  con¬ 
traction  has  been  applied  to  the  low  vibration  case. 
Also,  the  vortex  core  radius  was  redefined  to  be  equal 
to  the  vortex  core  measured  during  the  HART  pro¬ 
gram.  In  an  attempt  to  determine  the  eflfect  of  wake 
geometry  on  the  aeroacoustic  predictions  for  the  two 
cases  with  HHC,  the  orientation  of  the  wake  with  re¬ 
spect  to  the  tip  path  plane  was  adjusted  so  that  the 
miss  distance  between  the  blade  and  the  vortex  would 
be  similar  to  the  measurements  made  of  the  miss  dis¬ 
tance  during  the  HART  program.  In  addition  for  the 
low  vibration  case,  a  double  vortex  structure  was  used 
as  measured  during  the  HART  test.  The  predicted, 
acoustic-pressure  time  histories  for  the  three  test  cases 
at  the  advancing-side  microphone  locations  listed  in 
Table  2  are  shown  in  Fig.  4.  The  retreating-side  mi¬ 


crophones  are  shown  in  Fig.  5.  The  top  row  contains 
the  pre-test  predictions  and  the  bottom  row  contains 
the  post-test  predictions.  For  the  baseline  case,  the 
wake  re-structuring  improves  the  peak-to-peak  levels 
of  the  wave  form  on  the  advancing  side.  Not  only  is  the 
peak-to-peak  level  improved  on  the  advancing  side  for 
the  low  noise  case,  but  the  overall  wave  form  compares 
better  to  the  measured  data.  The  changes  to  the  wave 
forms  due  to  the  wake  re-structuring  do  not  improve 
the  predictions  for  the  low  vibration  case.  Perhaps  the 
wake  re-structureing  needs  to  be  optimized  for  each 
test  case.  The  contour  plots  of  the  mid- frequency  SPLs 
for  the  three  test  cases  are  shown  in  Fig.  6.  The  con¬ 
tour  plots  in  the  first  row  are  the  pre-HART  predic¬ 
tions  and  the  contour  plots  in  the  second  row  are  the 
post-HART  predictions.  The  wake  re-structuring  only 
slightly  affects  the  contour  plots  for  the  baseline  case. 
However,  it  improves  both  the  magnitude  and  the  di¬ 
rectivity  characteristics  of  the  low  noise  and  the  low 
vibration  cases. 

4.2  NASA  LaRC  Pre-HART  and  Post-HART 
Acoustic  Predictions 

The  methodology  used  to  predict  the  spanwise  blade 
loading  for  the  pre-HART  acoustic  predictions  made 
by  NASA  LaRC  is  also  discussed  in  detail  in  Ref. 
6.  The  prediction  method  begins  with  a  10  degree 
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Fig.  4  DLR  predicted,  acoustic-pressure  time  histories  for  the  advancing  side  microphone  locations  listed  in 
Table  2.  Row  1)  pre-test  predictions.  Row  2)  post-test  predictions. 
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Fig.  5  DLR  predicted,  acoustic-pressure  time  histories  for  the  retreating  side  microphone  locations  listed  in 
Table  2,  Row  1)  pre-test  predictions.  Row  2)  post- test  predictions. 
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Fig.  6  DLR  predicted,  mid-frequency,  SPL  contours.  Row  1)  pre-test  predictions.  Row  2)  post-test  predicitons. 


azimuthal  resolution  of  the  rotor  wake  calculated  by 
CAMRAD.Modl.  HIRES  then  computes  the  air¬ 
loads  at  1  degree  azimuthal  increments  and  at  75 
radial  stations.  The  rotor  blade  dynamics  are  calcu¬ 
lated  by  CAMRAD.Modl.  For  the  post-HART  pre¬ 
dictions,  the  measured  blade  dynamics  are  used  to 
prescribe  the  blade  motion  in  CAMRAD.Modl  before 
the  HIRES  computation  of  the  airloads.  The  NASA 
LaRC  predicted,  acoustic-pressure  time  histories  at 
the  advancing-side  microphone  locations  for  the  three 
test  cases  are  shown  in  Fig.  7.  The  retreating-side  mi¬ 
crophone  locations  are  shown  in  Fig  8.  The  time  histo¬ 
ries  in  the  top  row  are  the  pre-test  predictions  and  the 
time  histories  in  the  bottom  row  are  the  post  test  pre¬ 
dictions.  The  peak-to-peak  levels  are  reduced  on  the 
retreating  side  of  the  rotor  for  the  baseline  case  when 
the  measured  blade  dynamics  are  used  for  the  predic¬ 
tion,  but  the  waveform  does  not  change.  For  the  low 
noise  case,  the  low  frequency  character  of  the  waveform 
is  increased  and  the  peak-to-peak  levels  are  decreased 
by  the  application  of  the  measured  blade  dynamics. 
The  low  vibration  case  is  similar  to  the  baseline  case, 
wherein  the  peak  to  peak  levels  are  reduced  on  the 
retreating  side  without  much  change  to  the  waveform 
when  the  predicted  blade  dynamics  are  replaced  by 
the  measured  blade  dynamics  in  the  prediction  of  the 
airloads.  The  pre-test  and  post-test,  LaRC-predicted 
contour  plots  of  the  mid-frequency  SPLs  are  shown  in 
Fig.  9.  For  the  baseline  case,  the  SPLs  are  reduced 


in  the  post-test  calculations.  Also,  the  local  peak  on 
the  retreating  side  has  been  eliminated.  The  elimina¬ 
tion  of  the  local  peak  in  the  post-test  prediction  is  also 
evident  for  the  low  noise  case.  For  this  case,  the  lev¬ 
els  over  the  whole  contour  have  been  decreased.  The 
low  vibration  case  changes  only  slightly  by  the  appli¬ 
cation  of  the  measured  blade  motion  to  the  airloads 
prediction.  The  peak  levels  on  the  retreating  side  are 
reduced  without  the  elimination  of  the  retreating  side 
local  maximum. 

4.3  ONERA  Pre-HART  and  Post-HART 
Acoustic  Predictions 

A  detailed  description  of  the  pre-test  ONERA  predic¬ 
tion  methodology  can  be  found  in  Ref.  6.  The  pre¬ 
diction  methodology  consists  of  a  vortex  sheet  roll¬ 
up  model  that  results  in  interacting  vortices.^®  For 
the  pre-HART  predictions,  the  viscous  vortex  core  ra¬ 
dius  in  ARHIS  was  determined  by  an  evolution  law 
based  on  extrapolated  hot-wire  measurements  of  a 
two-bladed  rotor  model  in  hover.^  The  viscous  core 
radii  obtained  with  this  law  are  often  smaller  than  0.1 
chord  on  the  advancing  side  and  approximately  0.15 
on  the  retreating  side.  The  post-test  AHRIS  vortex 
core  size  model  uses  an  improved  evolution  law  based 
on  the  preliminary  analysis  of  the  HART  LDV  mea¬ 
surements.  The  adjusted  law  results  in  viscous  core 
radii  of  0.25  chords  on  the  retreating  side  that  are  in 
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Fig.  8  NASA  LaRC  predicted,  acoustic-pressure  time  histories  for  the  retreating  side  microphone  locations 
listed  in  Table  2.  Row  1)  pre-test  predictions.  Row  2)  post-test  predictions. 


Baseline  (Run  140)  Low  Noise  (Run  138)  Low  Vibration  (Run  133) 


Fig.  9  NASA  LaRC  predicted,  mid- frequency,  SPL  contours.  Row  1)  pre-test  predictions.  Row  2)  post- test 
predicitons. 


agreement  with  the  experimental  data.  The  viscous 
core  radii  on  the  advancing  side  only  changes  slightly. 
The  ONERA  predicted,  acoustic-pressure  time 
histories  of  the  three  test  cases  are  shown  in  Fig.  10  for 
the  advancing  side  microphone  locations.  The  retreat¬ 
ing  side  microphone  locations  are  shown  in  Fig.  11. 
The  only  difference  between  the  pre-HART  (top  row) 
and  the  post-HART  (bottom  row)  acoustic-pressure 
time  histories  for  aU  three  test  cases  is  the  magnitude 
of  the  peaks.  This  is  the  expected  result  of  increasing 
the  core  radius  in  the  airloads  computation.  The  ON¬ 
ERA  predicted,  pre-test  and  post-test,  mid-frequency 
SPLs  are  shown  in  Fig.  12.  The  only  difference  be¬ 
tween  the  pre-test  and  post-test  SPL  contours  are  the 
lower  levels  for  the  post-test  predictions. 

4.4  AFDD  Pre-HART  and  Post-HART 
Acoustic  Predictions 

The  main  differences  between  the  pre-HART  predic¬ 
tions  and  the  post-HART  predictions  made  by  AFDD 
are  a  change  to  the  partial  angle  calculation  and  the 
fine  tuning  of  the  dynamic  model  of  the  rotor  blade 
in  an  attempt  to  improve  the  dynamic  blade  motion. 
The  AFDD  pre-HART  prediction  methodology  is  sum¬ 
marized  in  Ref.  6.  A  closer  look  at  the  partial  angles 
calculated  by  C  AMR  AD/ J  A  revealed  that  the  angles 
were  not  being  calculated  correctly.  C  AMR  AD/ J  A 


computes  the  partial  angles  by  subtracting  the  induced 
velocity  caused  by  the  tip  vortices  within  the  compu¬ 
tational  domain  of  FPR  from  the  inflow  angles.  Spec¬ 
ifying  a  core  radius  of  2.0  chords  for  the  partial  an¬ 
gle  calculation  eliminates  the  anomalies  from  the  par¬ 
tial  angle  calculation.  A  core  radius  of  0.20  chords 
is  used  for  all  other  CAMRAD/JA  and  FPR  calcu¬ 
lations.  The  AFDD  predicted,  acoustic-pressure  time 
histories  of  the  three  test  cases  are  shown  in  Fig.  13 
for  the  advancing-side  microphone  locations.  The  pre¬ 
dicted  acoustic  pressure  for  the  retreating-side  micro* 
phone  locations  are  shown  in  Fig.  14.  There  is  only  a 
slight  change  in  magnitude  between  the  pre-test  and 
post-test  predictions  without  much  change  in  wave¬ 
form  for  both  the  advancing  and  retreating  side  micro¬ 
phone  locations.  The  AFDD  predicted,  mid-frequency, 
SPL  contours  are  shown  in  Fig.  15.  There  is  little  dif¬ 
ference  between  the  pre-test  and  post-test  prediction 
of  the  SPLs  for  the  baseline  case.  The  slight  decrease 
in  contour  levels  for  this  case  is  most  likely  caused  by 
the  change  in  the  partial  angles.  The  contour  levels 
for  the  low  noise  case  actually  increase  for  the  post- 
HART  predictions  with  only  a  slight  change  in  direc¬ 
tivity.  The  actual  contour  levels  for  the  low  vibration 
case  change  only  slightly,  but  the  local  maximum  on 
the  advancing  side  has  moved  aft  for  the  post-test  pre¬ 
dictions. 
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Fig.  10  ONER  A  predicted,  acoustic-pressure  time  histories  for  the  advancing  side  microphone  locations  listed 
in  Table  2.  Row  1)  pre-test  predictions.  Row  2)  post-test  predictions. 
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Fig.  11  ONERA  predicted,  acoustic-pressure  time  histories  for  the  retreating  side  microphone  locations  listed 
in  Table  2.  Row  1)  pre-test  predictions.  Row  2)  post-test  predictions. 
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Fig.  12  ONERA  predicted,  mid-frequency,  SPL  contours.  Row  1)  pre-test  predictions.  Row  2)  post-test 
predicitons. 
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Fig.  13  AFDD  predicted,  acoustic-pressure  time  histories  for  the  advancing  side  microphone  locations  listed  in 
Table  2.  Row  1)  pre-test  predictions.  Row  2)  post-test  predictions. 
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Normalized  time  Normalized  time  Normalized  time 

Fig.  14  AFDD  predicted,  acoustic-pressure  time  histories  for  the  retreating  side  microphone  locations  listed  in 
Table  2.  Row  1)  pre-test  predictions.  Row  2)  post-test  predictions. 


Baseline  (Run  140)  Low  Noise  (Run  138)  Low  Vibration  (Run  133) 


Fig.  15  AFDD  predicted,  mid- frequency,  SPL  contours.  Row  1)  pre-test  predictions.  Row  2)  post-test  predici- 
tons. 
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Fig.  16  Predicted,  acoustic-pressure  time  histories  using  a  common  set  of  aerodynamic  input  for  Rim  133.  Row 
1)  advancing  side  microphone  location  listed  in  Table  2  for  Run  133.  Row  2)  retreating  side  microphone  location 
listed  in  Table  2  for  Run  133.  Row  3)  in-plane  microphone  location. 

5  ACOUSTIC  PREDICTIONS  USING  determine  the  blade  attitude.  Also,  the  uncorrected 

COMMON  AERODYNAMIC  INPUT  shaft  angle  was  used  to  define  the  location  of  the  blade 

with  respect  to  the  observer.  The  resulting  predicted, 
acoustic-pressure  time  histories  are  shown  in  Fig.  16 
Even  though  all  the  prediction  team  members  use  a  nu-  advancing  and  retreating  side  microphone  lo- 

merical  simulation  of  the  Ffowcs  WiUiams  and  Hawk-  ^  2  for  Run  133  and  for  a  microphone 

ings  equation,  each  implementation  is  different.  To  ^his  in-plane  micrc^ 

determine  if  there  are  any  differences  in  the  acous-  mounted  on  the  nozzle  and  is  at  x  =  -6.71  m, 

tic  predictions  caus^  by  the  specific  implementations  ^  ^  3  395  m,  and  z  =  0.0.  There  are  only  a  few  slight 

of  the  FW-H  equation,  all  the  prediction  team  mem-  differences  in  magnitude  for  the  acoustic  predictions 

bers  used  a  common  set  of  aerodynamic  blade  surface  advancing  side  microphone.  The  differences  in 

pressures.  The  common  set  is  a  prediction  of  the  low  ^j^tude  are  much  greater  for  the  predictions  at  the 

vibration  case  (Run  133)  provided  by  ONERA.  Each  retreating  side  microphone.  For  the  predictions  at  the 

team  member  used  the  blade  loading  at  the  same  30  microphone,  the  main  difference  is  the  magni- 

radial  stations  and  41  chordwise  positions  for  the  non-  evident  in 

compact  model  of  the  dipole  source.  All  acoustic  pre-  predictions.  The  predicted,  mid-frequency, 

dictions  for  this  case  were  made  assummg  that  the  gp^  contours  are  shown  in  Fig.  17.  The  character 
rotor  was  stiff,  i.e.  no  flapping,  bendmg  or  torsion.  directivity  pattern  is  quite  similar  for  all  four 

The  measured  cyclic  and  collective  pitch  were  used  to 
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Fig.  17  Predicted,  mid-frequency,  SPL  contours  for  Run  133,  using  a  common  set  of  aerodynamic  input. 


methods.  NASA  LaRC  predicts  higher  levels  down¬ 
stream  of  the  rotor  than  do  the  rest  of  the  prediction 
team,  and  the  local  maximum  on  the  advancing  side 
is  slightly  higher  for  the  ONERA  predictions. 

6  SUMMARY  AND  CONCLUSIONS 

The  pre-test  and  post-test  acoustic  predictions  for 
the  Higher  Harmonic  Control  Aeroacoustic  Rotor  Test 
have  been  presented.  The  comparisons  between  the 
pre-test  and  post-test  predictions  indicate  how  the 
HART  program  can  guide  researchers  to  improve  the 
aeroacoustic  and  dynamic  modeling  of  the  BVI  phe¬ 
nomena  so  that  low  noise  concepts  such  as  HHC  can 
be  explored.  To  improve  their  prediction  capability, 
the  DLR  combined  a  contraction  of  their  prescribed 
wake  model  with  a  re-orientation  of  the  wake  geometry 
with  respect  to  the  tip  path  plane.  This  re-structuring 
of  the  wake  can  improve  the  aeroacoustic  predictions 
and  will  be  used  to  refine  the  DLR  free- wake  model. 
The  substitution  of  the  measured  blade  dynamics  for 
the  predicted  dynamics  does  improve  the  aeroacous¬ 
tic  predictions  made  by  NASA  LaRC.  This  example 
demonstrates  the  necessity  of  properly  calculating  the 
blade  motion  for  BVI  noise  predictions,  especially  for 
cases  with  HHC  input.  ONERA ’s  post-test  predic¬ 
tions  show  that  changing  the  core  size  of  the  interact¬ 
ing  vortices  to  more  closely  agree  with  the  measured 
core  sizes  improves  the  prediction  of  the  magnitude 
of  the  predicted,  acoustic-pressure  time  histories  and 
the  SPL  contours.  AFDD’s  fine  tuning  of  the  dynamic 
model  of  the  rotor  blade  affected  the  aeroacoustic  pre¬ 
dictions  only  slightly.  The  prediction  of  the  dynamic 
blade  motion  by  CAMRAD/JA  is  still  not  adequate 
for  the  cases  with  HHC  input. 

There  are  distinct  diflferences  in  the  acoustic  predic¬ 
tions  made  by  each  of  the  prediction  team  members 


when  the  same  aerodynamic  input  is  used.  The  exact 
cause  of  these  difference  is  not  currently  obvious,  but 
documenting  the  extent  of  the  differences  should  help 
the  rotorcraft  community  imderstand  where  the  state- 
of-the-art  lies  in  acoustic  prediction  methodology. 
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